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Abstract. For MSSM phenomenology, soft SUSY breaking dimension-three operators are impor- 
tant, in particular the couplings between Higgs bosons and squarks. In scattering processes, per- 
turbative unitarity is violated at modest center-of-mass energy if these couplings are much larger 
than the masses of the scalar particles involved. Assuming perturbative unitarity, constraints on 
the trilinear couplings can be determined using a computer program that we have developed. 
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1 Introduction 

For phenomenology in the Minimal Supersymmetric 
Standard Model (MSSM), the soft supersymmetry 
(SUSY) breaking dimension-three trilinear couplings 
between Higgs bosons and squarks play an important 
role. They determine the mass splitting in the stop 
and sbottom sector and can have a profound influence 
on Higgs physics. Unfortunately, these SUSY breaking 
couplings are essentially free parameters. Thus, any 
constraints that can be placed by theoretical consider- 
ations are useful. Such constraints can be obtained, for 
example, by considering the symmetry breaking of the 
scalar sector and requiring that there are no color or 
charge breaking minima pTj. Requiring that the MSSM 
allows for a perturbation theory treatment, we here 
consider a different approach: perturbative unitarity 
in 2 — > 2 scattering of scalar particles limits the size of 
the trilinear couplings. 

The approach is similar to the consideration of lon- 
gitudinal weak boson scattering in the Standard Model 
(SM), where the absence of perturbative unitarity vi- 
olation at high energies provides an upper bound on 
the Higgs quartic coupling and thus on the Higgs boson 
mass [2l|3] . For trilinear couplings in the MSSM scalar 
sector, these unitarity violations arise at intermediate 
energies, somewhat above pair production thresholds, 
if trilinear couplings are chosen much larger than the 
masses of the scalar particles in the scattering process. 
The approach is applied to the trilinear couplings be- 
tween Higgs bosons and 3rd generation squarks, given 

by 

which is part of the MSSM soft SUSY breaking poten- 
tial. Here Ql is the 5[/(2)i-doublet of 3rd generation 
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squarks, b]^ and t]^ denote the right-handed sbottom 
and stop singlet fields, and Hd and Hu are the Higgs 
boson doublet fields. We want to limit the dimension- 
ful parameters which multiply the Yukawa couplings 
Ab = mfj/w cos and Xt —mt/v sin /S, and fi, the Higgs 
mixing parameter. 

2 Perturbative unitarity for scalars 

The starting point is the unitarity of the 5* matrix, 

S-tfi* = 1 , or cquivalently - i {T - T^) = T'^T , 

for the transition operator T with 5'= l + iT. To eval- 
uate this equation one usually restricts to 2 — > 2 scat- 
tering (which is a good approximation in perturbation 
theory) and then uses angular-momentum conserva- 
tion and symmetries of the model to partially diago- 
nalize T . In our analysis we additionally restrict our- 
selves to scalar fields. Let {f\T\i) denote the transition 
matrix elements with initial state \i) and final state 
I/), and let Tfi be the matrix element obtained from 
Feynman diagrams in momentum representation, 

{2nf 5(4) (p^ „ p^) t^^(^^ eos0) = {f\T\i) , 

evaluated in the center-of-mass system, where \/s is 
the total energy and 9 is the scattering angle. This 
transition is defined by the (ordered) particle content 
of the two states, as well as and 9, and depends 
on the masses m;^ of the particles in state l — i,f via 
the functions A; = A(s, to^-^, 771^2); X{x,y, z) — x'^ + y'^ + 
z'^ — 2xy—2yz~2zx. The dependence on the scattering 
angle 9 is eliminated by projection onto partial waves 
of total angular momentum J =0, 1, 2, . . . 

^fi = o ' / dcos 9 ff, (Ti, cos 9) Pj (cos 9) . 
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Fig. 1. Scalar 2— >2 tree level scattering diagrams. 

The Pj are the Legendre polynomials. The factor 1/2 
is a standard convention and leads to the factor 1/2 in 
Eq. 11]). In this normalization, extra factors of l/\/2 
have to be included for each state with two identi- 
cal particles. Higher partial waves usually give smaller 
amplitudes, so only J = 0,1 amplitudes have to be con- 
sidered in practice. The unitarity condition now reads 

^ {T/^~T^)=Y.'^i!;Ti!^. (1) 
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The sum is taken over intermediate states. Restriction 
to only relevant and two-particle scalar states in the 
sum slightly underestimates the right-hand side and 
leads to conservative bounds. 

The 'true', physical matrix is normal and can 
therefore be diagonalized. The same holds for the Born 
amplitude, which we use instead. The diagonalized ma- 



trix and thus the eigenvalues satisfy 



ImT/ ^ IT 



J|2 



(2) 



The 'true' eigenvalues (for any given energy ^/s) must 
lie on the circle (called Argand diagram) given by ([2|) : 
y = x^+y^ for x =Re'7^f and y =Im7j/ which implies 
\x\ < 1/2. For the Born approximation, the phases of 
the fields can be chosen such that all 2^2 amplitudes 
are (nearly) real, if CP (nearly) holds. Approximat- 
ing X by the corresponding Born amplitude yields the 
desired unitarity bound. The circle restricts |a;| <l/2, 
which is the unitarity bound generally used in the lit- 
erature. A Born value ofa;= 1/2, y = needs at leaslQ 
a correction of -\/2— 1~41% to become unitary. Such 
large corrections indicate a breakdown of perturbation 
theory. In addition to the perturbative unitarity bound 
of I a; I < 1/2 we therefore also consider |a;| < 1/6 as a 
condition for which the Born amplitude remains sufh- 
ciently small to trust perturbation theory. 



2.1 A toy model 

Figure [T] shows generic tree level Feynman graphs for 
ipiip2 — *■ </'3V'4 scattering with trilinear couplings of 
(possibly different) scalar fields ipi (/ = 1,.., 5). The cor- 
responding amplitudes are A^/(q^— '^s)? where = 
s, i, M is one of the Mandelstam variables, and A stands 
for the trilinear couplings, t and u depend on ^/s, cos 9, 
and the masses to/ of the exterior particles. Project- 
ing onto partial waves, the amplitude for J = has a 
structure roughly like: 



j=o 



fi. 
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max{s, to|} ' 



(3) 



^ 'At least' means that this would be the minimal case 
where the correction directly hits the nearest circle point. 



where a factor of 2 has been assumed to account for the 
partial wave projection. The second factor is smaller 
than 1 and the third factor becomes large if both the 
scattering energy ^/s and the mass of the internal par- 
ticle are small compared to the couplings A. Highest 
values are usually found for energies near the kine- 
matic threshold, i.e. at energies where the model should 
work properly, in contrast to weak boson scattering in 
the SM 



2.2 Handling poles 

Clearly the Born amplitudes are not sufficient to de- 
scribe scattering processes where intermediate parti- 
cles become on-shell. In the s channel in Figure [T] this 
happens when ^/s = tos . Born amplitudes only will be 
used for unitarity considerations and s-channel poles 
are cut out by the condition 



> am r(Q = bm) 



(4) 



Here a,b>l are (suitably chosen) constants and the 
'running width' r{Q) of the internal particle ip at en- 
ergy Q is approximated by the decay width via replac- 
ing its mass to with the energy Q in the phase space 
factor. This condition ^ has to be fulfilled for all in- 
ternal particles appearing in the s channel. If this is 
not the case, the amplitude is set to zero, as well as 
the irreducible part of Tj[^'^ this process is in, because 
of possible destructive interference of matrix elements. 

A width cannot be included in the Born propagator 
for two reasons: First, Tfi is no longer diagonalizable 
(at this level of approximation). Second, our r{^/s) 
grows (linearly for large y^) with ^/s, which is not a 
good approximation to the propagator as y/s^m. 

The internal particle in the u channel of Figure 
[1] can also become on-shell for certain combinations 
of masses. This occurs e.g. if tpi can decay into (^4, tp^ 
and ip2, (p5 can fuse to ips (MSSM example: t2ti — > t2ti 
with u channel /i" exchange when m^^ > TO^^-f to/j). One 
obtains another possibility by switching labels 1 <-> 2 
and 3 4 or two similar conditions for the t channel 
by exchange of 3^4. The first case has the condition: 

CTOi > m4 + TOs A m2 + TO5 < CTO3 . (5) 

with some suitably chosen constant c > 1. Amplitudes 
where a condition like in ([5]) is fulfilled cannot be com- 
puted because the internal particle becomes on-shell 
for some value of the scattering angle. The constants 
a, 5, c are chosen larger than one because in proximity 
of a pole one encounters unphysical enhancements of 
the amplitude in the pure Born approximation. Still, 
some enhancement can appear in special cases. 

If some Born matrix elements cannot be calculated 
because of a t or u channel pole, the tree level matrix 
Tj[ cannot be diagonalized. The solution is a partial 
diagonalization. Assuming time reflection invariance, 
we write the left-hand side of ^ as ImTy^. Define the 
set B of all kinematically accessible states at given ^/s 
from an irreducible part of 7y{ and the set CcB such 
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that h,lGC satisfy at or u channel pole condition for 
h^l. Then, for states /, i€B\C equation ([1]) becomes 

h<£B\C h(£C 

We diagonahze the sub-matrix ('?}i)|/,ieB\c of all 
states accessible at a given energy ^/s with a unitary 
matrix U and write the diagonalized part as T to ob- 
tain: 

hec 

The second term on the right-hand side (denoted by 
R^) is positive. At low energies where only a few states 
are kinematically allowed, typically vanishes. The 
circle equation y — x^+ip now holds for y — lmT^f and 

\x\ = {{Refdf+R^Y'\ (6) 
3 MSSM constraints 

The MSSM has a lot of (> 100) parameters, many of 
which are constrained by measurement of the prop- 
erties of SM particles and (nearly) conserved symme- 
tries. Here we focus on some of the SUSY breaking 
parameters, namely the interplay of mass terms for 
the scalar fields and the dimensionful A parameters 
of the trilinear couplings, and also /i. These param- 
eters are constrained by lower mass bounds for the 
sparticles and the Higgs bosons [S]. Their phases are 
bounded by the absence of large CP violation. In addi- 
tion, the fields must be in the observed minima of the 
potential and have positive masses after electroweak 
symmetry breaking. Complex, strong constraints from 
not being in charge and color breaking minima were 
derived in [Tj. The new constraints from perturbative 
unitarity should be considered in addition. 

3.1 Particles, parameters and calculation 

In the MSSM large trilinear couplings may appear be- 
tween 3rd generation squarks and the Higgs bosons 
and the longitudinal degrees of freedom of and 
To work in the scalar sector only, the Goldstone 
equivalence theorem and the Feynman i?^=i-gauge are 
used, such that the Goldstone bosons G^, G° represent 
the longitudinal polarizations of W^^ . The trilinear 
coupling strengths are given by vertex factors like 

V{G°, II, iji) = g' {At - ti* cot P) mt/2mw ■ 

Other stop vertices with Higgs bosons are similar and 
contain combinations of At and /i (in decoupling sce- 
narios mainly Xt= At — n* cot f3) . Likewise two sbot- 
tom Higgs boson couplings grow with At and fi, but 
are suppressed compared to the stop case, because 
they contain the factor mh/2mw- Some mixed stop- 
sbottom couplings also behave like the pure stop case, 
e.g. the vertex V{G~^,bL,iji) = i g' Xt mt/V^mw- 



Relevant for our analysis are the scalar Higgs bos- 
ons ft." H^, A^, H^, G°, and G^ and the heavy squark 
mass eigenstates bi, &2, ii, and i2- In addition to SM 
parameters, results will depend on the squark mass 
parameters M"^, M~, and M~ , the mass of the pseudo- 
scalar Higgs boson ttia, the Higgs mixing parameter fj,, 
the ratio of Higgs bosons VEVs tan (3 and the trilinear 
coupling parameters Ab and At- Neglecting any CP- 
violating effects, these parameters are chosen real in 
the following. 

With the scalar fields listed above, two-particle 
states are formed. One uses charge, color, and the 
non-existence of three-squark-vertices to form 15 inde- 
pendent blocks in the scattering matrix. The biggest 
block (charge 0, color singlet) contains 21 states, at 
high enough energy, and usually supplies the largest 
eigenvalue. A Monte Carlo search over the scattering 
energy is used to find the largest eigenvalue of the 
scattering matrix, which is obtained by numerically di- 
agonalization. We use, by default, MS quark masses (at 
120 GeV) to calculate the squark vertices and masses. 

Large trilinear couplings can lead to strong loop 
suppressions of SUSY Higgs masses 6^ in regions where 
our analysis indicates a breakdown of perturbation 
theory. In such cases we ignore clashes with LEP 
bounds [5] and conservatively use m^o = 120 GeV in- 
stead. 



3.2 Relevance of the method for the MSSM 

Because only heavy squarks and Higgs bosons were 
used, no parameters for charginos and other squarks 
or sleptons were needed for the calculation. The results 
are also practically independent of these other param- 
eters. The maximum mixing case, > 2 times the 
sfermion mass scale, is often used as a benchmark sce- 
nario. An important bound for these cases often used 
in literature derives from [4| (see e.g. Eq.(5) in [1]): 

A? < 3 ( M| + M? + -f l/ip) , (7a) 
A?<3(M?+M2 + ml,^+|/.|2), (7b) 

with rr?^ = cos^/3 (sin^/3) ± 1/2 to^ cos 2/3. 

For large tan/3 this means |Atp <3 (M?-|-ikr?-m|/2) 
and l^bp < 3 {MyrW^^m^j^rn^^l^). The At bound in 
([7a|) is mostly determined by M% + M~. The unitarity 
bounds instead depend on the masses of the lightest 
particles with large trilinear couplings, which there- 
fore can be stronger in the case of exceptionally light 
scalars while retaining a large value for Af^ -I- M?. 

3.3 Examples for unitarity bounds 

The approach of the unitarity bound for large trilin- 
ear couplings is illustrated in the figures. MSSM pa- 
rameters for these two examples are given in the cap- 
tions. Figure [2] strongly constraints /z, because A^ and 
At have been set only slightly below the bound from 
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Fig. 2. ^ bound example with tan/3 = 50, mA = 100 GeV, 
Mq=400 GeV, Mi = 600 GeV, Mj = 630 GeV, ^t = l TeV, 
and ^6 = 1.3 TeV. 135 GeV < mj^ < 404 GeV varies witli /i. 
\x\ is the maximal eigenvalue of T^^ as modified in Eq. (6). 
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Fig. 3. A bound example, together with the bound from 
Eq. (|7a|l (vertical lines) using tan/3 = 30, = 1 TeV, 
M = -200 GeV, Mq=Mi = 2.9 TeV with fixed = 100 
GeV and m^o = 120 GeV. 



Eq. ([7]) and a high tan/3 has been chosen. Alterna- 
tively, one could consider Ab bounds for large, fixed 
H in this scenario. In Figure [2] the diagonalized am- 
plitude mainly results from the interfering low energy 
processes bibl ^ h°h°, ^ h°H°, bill ^ bibl, and 
bibl^H°H°. 

Figure [3] shows an example for the case of a light 
ii held at a fixed mass of 100 GeV by adjusting M? 
for each value of ji. To evade the bounds ([7]) for A = 
Ah — At, Mq is chosen high, which leads to the max- 
imum mixing case, where the ti,t\ couplings to the 
Higgs bosons are strongly reduced. The largest am- 
plitude is almost completely due to the tit\ h^h^ 
i-channel ti exchange diagram and is mostly indepen- 
dent of tan/3 ^ 1, niA 3> m^, ^ Mq and M^. In 
this scenario Eq. ([7a|) gives |^| ^ 5 TeV, while the 
perturbativity condition < 1/6 yields |A| <4.4 TeV. 

4 Conclusion 

A general method has been developed for constraining 
trilinear couplings of scalars by using perturbative uni- 
tarity. It has been implemented for the MSSM, more 
specifically for third generation squarks and Higgs bos- 
ons, but it can also be used for other models. Since 
trilinear couplings of scalars correspond to superrenor- 
malizable dimension three operators, they produce the 
largest contributions to scattering amplitudes at mod- 
est center-of-mass energy. Thus, perturbative unitarity 
bounds are derived from scattering amplitudes at ener- 
gies not very far above production threshold and they 
are therefore independent of the ultraviolet structure 
of the theory. One finds that trilinear couplings cannot 
be much larger than the masses of the scalar particles 
involved. Quantitatively, ratios of more than about a 
factor 5 are forbidden (see Eq. ([3]) ) . 

We have analyzed two MSSM scenarios with small 
sbottom or stop masses, trying to constrain the soft 
SUSY breaking parameters A^ and At and the Higgs 
mixing parameter /i at large tan/3. Mixing effects be- 
tween left- and right-handed squarks lead to a com- 



plicated picture, since large A parameters may lead to 
only modest trilinear couplings of the lightest scalars 
when the relevant mixing angle factors are small. We 
have built a numerical program to study these effects 
for two Higgs doublets and third generation squarks 
within the MSSM for arbitrary input parameters tan /3, 
/i, ruA, A}), At, M^, M|, and M?. For small ^, result- 
ing constraints are similar in strength to bounds de- 
rived from the exclusion of false vacua e.g. Eq. ([7]), 
however, the underlying arguments are quite differ- 
ent. Our approach therefore provides a complementary 
method for constraining soft SUSY breaking param- 
eters. Given a specific MSSM scenario, specified by 
the Lagrangian parameters, our program gives a test 
whether these parameters clash with unitarity. 

Since supersymmetric theories are at heart pertur- 
bative, even stronger constraints are obtained by the 
requirement that scattering amplitudes sufficiently far 
from resonances do not leave the perturbative domain. 
Such a test is also possible with our program, e.g. by 
requiring that the largest eigenvalue of the scattering 
matrix stays well below the unitarity bound of 1/2. 
Such an analysis becomes instructive when analyzing 
higher loop effects involving SUSY particles which po- 
tentially lead to large corrections of parameters for SM 
particles, like Higgs-fermion Yukawa couplings or the 
mass of the lightest CP-even Higgs boson. We have 
not yet performed such a systematic analysis. 
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